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Abstract. We investigate the parametric fluctuations in the quantum survival 
probabiUty of an open version of the ^-kicked rotor model in the deep quantum 
regime. Spectral arguments [Guarneri I and Terraneo M 2001 Phys. Rev. E 65 
01 5203 (R)] predict the existence of parametric fractal fluctuations owing to the strong 
dynamical localisation of the eigenstates of the kicked rotor. We discuss the possibility 
of observing such dynamically-induced fractality in the quantum survival probability 
as a function of the kicking period for the atom-optics realisation of the kicked rotor. 
The influence of the atoms' initial momentum distribution is studied as well as the 
dependence of the expected fractal dimension on finite-size effects of the experiment, 
such as finite detection windows and short measurement times. Our results show 
that clear signatures of fractality could be observed in experiments with cold atoms 
subjected to periodically flashed optical lattices, which offer an excellent control on 
interaction times and the initial atomic ensemble. 



PACS numbers: 05.45.Mt, 42.50.Vk, 05.60.Gg 
1. Introduction 

Experiments witli cold atoms nowadays offer unique possibilities for the study of 
single particle motion and collective particle dynamics in tailored optical or magnetic 
potentials. The atomic centre-of-mass motion can be prepared and controlled with 
unprecedented precision, what allows experimentalists to realise and study many toy 
models of condensed matter physics PJ. Since in experiments with cold atomic gases 
noise and perturbations can be driven to a minimum, which often is indeed negligible, 
such setups offer a great advantage with respect to solid-state realisations. 

In this paper, we discuss the possibility of observing sensitive quantum effects 
which manifest in a fractal variation of a transport function with respect to a well- 
tunable control parameter. Similar fractal fluctuations of the transmission probability 
across solid-state samples have been measured recently j2] in systems whose underlying 
classical phase space typically contains mixed regular-chaotic structures. Most features 
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of these experiments can be understood semiclassically as a consequence of the 
phase space topology El El- However, the precise origin of the observed fractal 
conductance fluctuations in these experiments is not yet fully understood [Ij, and, 
in fact, various theoretical models [3 |H] predict fractal conductance fluctuations for 
mesoscopic devices. Our aim is to design a concrete experimental scenario in which 
parametric fractal fluctuations could be measured with high precision cold-atom setups. 
In such experiments the cross-over between mixed and completely chaotic classical 
dynamics can be scanned easily jHl EH E] , and hence fractal transmission probabilities 
could be measured in a regime where classical or semiclassical arguments do not apply. 

As was shown by Guarneri and Terraneo [Zj , fractal fluctuations in the transmission 
probability of a quantum scattering problem arise naturally as a consequence of the 
spectral properties of the system. The two essential conditions on the spectrum are (i) 
a power-law distribution of decay widths and (ii) uncorrelated real parts of the energy 
spectrum. Moreover, various eigenstates have to contribute together to the decay, a 
fact which is expressed formally by requiring that (iii) the average decay width is much 
larger than the mean level spacing. Based on these conditions, the theory of [7j explained 
the occurrence of quantum fractal fluctuations in the ^-kicked rotor model in the deep 
quantum realm |T2], where semiclassical arguments cannot explain the occurrence of 
fractality. 

In this paper we study a similar dynamical situation as in ^2] , yet with important 
modifications which fully account for the actual experimental realisation of the kicked 
rotor. Using either cold or ultracold atomic gases, the kicked rotor is realised by 
preparing a cloud of atoms with a small spread of initial momenta, which is then 
subjected to a one-dimensional optical lattice potential, fiashed periodically in time [T8] . 
Let us call ki the wave number of the optical lattice, f the fiashing period ("kicking" 
period), p the momentum of the single atom, x its centre-of-mass position, Vq the 
maximum potential depth, and M is the atomic mass. It is convenient to adopt rescaled 
units by noting that Pr = hk^ is the photon recoil momentum and Eji = {hkiY 1 2M is 
the recoil energy [01 HH Ej. So we define p = p / 2pR, x = x / 2kL, t = f ■ 8Er / h. The 
kicking strength of the lattice is expressed hj k = Vqt / h. The Hamiltonian now reads 
in dimensionless units [T5] 



Owing to the 5-interaction of the potential with the atoms, the time evolution operator 
between kicks can be explicitly written in a factorised form, extremely convenient for 
numeric simulations. The derivation of the one-period evolution operator exploits 
the spatial periodicity of the potential by Bloch's theorem [HI HH]. This defines 
quasimomentum /3 as a constant of the motion, the value of which is the fractional 
part of the physical momentum p in dimensionless units p = n + (3 (ra G N). Since 13 
is a conserved quantum number, p can be labelled using its integer part n only. The 
spatial coordinate is then substituted hy 9 = x mod(27r) and the momentum operator 
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In close analogy to the transport problem across a solid-state sample, we follow 
|12j to define the quantum survival probability as the fraction of the atomic ensemble 
which stays within a specified region of momenta while applying absorbing boundary 
conditions at the "sample" edges. If we call ip{n) the wave function in momentum space 
and ni < 712 the edges of the system, absorbing boundary conditions are implemented by 
the prescription ipiji) = if < ni or n > n2. This truncation is carried out after each 
kick. This procedure mimics the escape of atoms out of the spatial region where the 
dynamics induced by the Hamiltonian ((H) takes place. If we denote by P the projection 
operator on the interval ]ni,n2[ the survival probability after t kicks is: 



We will show in the following that signatures of fractality in the survival probability 
could be observed in modern atom-optical experiments, where the initial atomic 
ensemble has a finite, non-zero width in momentum space. In contrast to the work of 
[T^ , where the initial quasimomentum is scanned to arrive at the parametric observable 
-Psurv(/5), we investigate the behaviour of Psnwij) as a function of the best controllable 
parameter in the experiment, namely the time r which elapses between two successive 



After a brief review of the results of Guarneri and Terraneo |7j applied to the 
dynamically localised kicked rotor (section 12)), we discuss in section El our choices of 
the system parameters, which are guided by the experimental possibilities as well as 
the conditions stated in Our central results on the occurrence of fractal survival 
probabilities are presented for the limit of long-interaction times (section^]) as well as for 
experimentally accessible initial momentum distribution and interaction times (section 
EJ. Section ini finally concludes the paper. 

2. Conditions for fractal fluctuations of the survival probability 

Without a priori assumptions on the integrability or chaoticity properties of the classical 
analogue of the quantum system of interest, Guarneri and Terraneo jZj showed that 
fractal conductance fluctuations occur if certain conditions on the quantum spectrum 
of the open system are fulfilled. 

The first condition, (i) a power-law distribution of the decay widths, is indeed 
present in the weakly opened quantum kicked rotor ^H]- We verified this by 
diagonalising the one- kick evolution operator Up, after representing it in the basis of 
momentum states. The matrix was cut at the positions Ui and n2 to mimic the required 
absorbing boundary conditions. 

If either of the two cutoffs {ui or ^2) is chosen sufficiently large, the shape of 
the wavefunction in momentum space supports an exponential tail, independent of the 
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evolution time (after a short transit time oc k"^ at which dynamical localisation has fully 
developed [201 121] )■ For such a situation in the locahsed regime, the probability density 
of decay widths was found to be p(T) oc over more than 10 orders of magnitude 
in r, consistent with previous studies jHl dl 1221 123 121j • If, on the other hand, ni and 
n2 were decreased, dynamical localisation is gradually destroyed and the distribution 
deforms continuously, giving more weight to larger widths and less to the very small 
ones. Such a deformation was observed in the analogous context of ionisation rates of 
microwave-driven hydrogen Rydberg atoms j21]. Our choice of ni and n2 represents a 
compromise between the maximum width of typical experimental detection windows in 
momentum space and a guaranteed dynamically localised momentum distribution over 
a substantial interval of momenta. In the next section, we state the precise values of ni 
and n2 which we investigated in this paper. 

In the regime of strong dynamical localisation, the quasienergy spectrum of the 
5-kicked rotor has a Poisson-like statistics [23]. Under the same conditions as stated 
above on the cutoff values ui and n2, this property of the real-parts of the quasienergy 
spectrum remains even when the system is opened 1^ . Hence, also the second 
requirement for fractality of [Z], that (ii) the energy spectrum consists of uncorrelated 
sequences, is fulfilled in good approximation for the opened (5-kicked rotor in the presence 
of dynamical localisation. 

The third condition stated in [Zj is that the opening of the system is weak, but 
still sufficient to guarantee that (iii) the average decay width is much larger than the 
mean level spacing. For our choice of parameters and cutoff values ui and n2, also this 
condition of overlapping "resonance peaks" is fulfilled, as we verified numerically from 
the quasienergy spectrum of the truncated matrix representation of W/j. 

As exercised in [7], the conditions (i-iii) are sufficient to guarantee self-affine 
fiuctuations in the quantum survival probability, with a predicted fractal dimension D{ 
which is related to the exponent of the width distribution p(r) oc F^" by the following 
general formula D{ = 1 + a / 2 ^ 1.5 for a ^ 1. 

We repeat that parametric fractal fiuctuations in the survival probability of 
dynamically localised kicked rotor have already been found in J2], before their origin 
could be explained in [Z]. In this work, however, we scan a different parameter than 
the one used in jT2I, which corresponded to quasimomentum. Here we use the kicking 
period r as control parameter, which can be much better controlled in state-of-the-art 
experiments [El [HI El |IH| than the initial value of momentum HH UHl [23 [2S1 123 ■ On 
the other hand, the use of r confronts us with a new problem which is discussed in the 
following section. 
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3. Choice of parameters 

3.1. Dynamical localisation and classical chaos 

For our analysis the value of the kicking strength k was chosen in the range 2 ... 6, 
OT kr = 2. 8... 8. 4, going along with the transition from local to global chaos with 
increasing k in this range [201 1^- For our choice of kicking periods r = ^efr > 1 
jSni El], classical trajectories wandering about hierarchical structures of the classical 
phase space will not have a quantum analogue because those structures are too small 
to be resolved by the wave function. This means that the observed fluctuations indeed 
arise from quantum localisation effects and not from a semiclassical diffusion process. 



3.2. The kicking period as control parameter 

As reviewed in section |2l the sufficient conditions for the occurrence of fractal 
ffuctuations are fulfllled for choices of r for which the 5-kicked rotor exhibits dynamically 
localised behaviour. However, besides dynamical localisation the quantum 5-kicked rotor 
supports "quantum resonant" motion for speciflc values of r and quasimomentum (5 
EOj . Our goal is to avoid as much as possible the impact of the quantum resonances 
on the dynamics, such that we can clearly identify the origin of the fractality of the 
survival probabilities. Since the parameter we scanned is the kicking period r, we 
verifled that no signatures of quantum resonances are found in the analysed small range 
of T and for the applied, flnite kick numbers. 

The quantum kicked rotor shows ballistic growth of momentum, shortly a quantum 
resonance, if 

r G {47rs/g; s, g G N}, /? G {m/2s, < m < s; m, s G N}, (4) 

and in these cases the time dependence of energy on the number of kicks is [23 IHO] 

E{t; r) = rit^ + 0{t), with r/ ~ (fc/g)^^. (5) 

The denominator q in the rational factor of r is called the order of the resonance. The 
set containing all the resonances has zero Lebesgue measure in any interval of kicking 
periods, but we do care about it because the dependence on r of the survival probability 
Psurv(T; t) is continuous for a fixed, finite number of kicks and the fiuctuations we want 
to observe should be caused by dynamical localisation and not by quantum resonances. 
The preceding growth estimate (0) establishes that a resonance is suppressed for a time 
that increases more than exponentially with its order. One way to avoid contributions 
from the resonances is to use a judicious choice of the range of r and sampling grid G 
used for numerical simulations or experiments. We chose 

G = {Ti = To + i ■ 6T,i e {0,. . . ,m- 1}}, with m = 10\ (6) 

where the value of tq / 47r is a fraction of the golden mean: 

1^ = il£(v/5 - 1), s = 6142, q = 95403, ^ < 10"^^ . (7) 

47r lOg^ ^' ' ^ ' g 47r(v^-l) ^ ^ 
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For 6t = 9.98 x 10"'', also all other grid points in G are incommensurable to Air up 
to the used significant digits. We verified that the lowest order resonance in the range 
[Tmin, Tmax] ~ [1-4, 1.41] has q = 107 and that there is no crowding of resonances of 
order q < 2000 anywhere in this interval. Since we are not going to use times longer 
than 10^ kicks in our simulations, and kicking strength of order unity, the quadratic 
term is suppressed dramatically by the coefficient // in (jSJ, for all occurring resonances 
q > 107. Finally, we explicitly checked throughout the simulations that localisation is 
at work by inspecting the average energy and, for selected values of r, the shape of the 
wave function in momentum space, which shows a characteristic exponential decrease 
as explained below in section EiSl 

We also tried a quantitative approach for the choice of the grid along the r 
axis. If some resonance were important, any numerical selection method could detect 
it and prefer grids with points away from the quantum resonances. Our method is 
based on the maximisation in the "grids space" of a function F{G{tq,St)) that adds 
a contribution from each resonance, up to a maximum order, within a given interval, 
and this contribution is the larger the farther the resonance position in r [see Eq. (0}] 
is from the nearest point of the grid. This means that a "higher mark" is achieved by 
the grids whose points are away from the resonances. Formally we defined 

S = {Ans/q} n{q< g^ax} H [ 

''~miru ''~maxj 

F{GiTo,5T)) = J2 fr (min{ | - | ; G G} ) , 

(Ar) = At; fl'\Ar) = Ar'; f^'\Ar) = Ar/q,. . 

Different definitions of the weight function /,.(Ar) allow us to give more weight to 
resonances with smaller q > 107 (i.e., to those which infiuence the time evolution of 
a wider neighbourhood along the r axis). Of course, this programme requires detailed 
knowledge of the dynamics near the high-order resonances of g > 107, but this goal has 
not been theoretically accomplished yet. None of our weight functions could resolve the 
presence of a resonance by a sharp minimum when applied to a specific grid. 

As a consequence of our choice of the interval of kicking periods and grid points 
in this interval, no signatures of quantum resonances are expected to manifest for 
interaction times of up to 10^ kicks. 

3.3. The opening of the system 

The probability decay arises from the open geometry of our system, which is 
implemented mathematically by imposing absorbing boundary conditions in momentum 
space [12]. This means that 

ilj[n) =0 if n < Ui < or n > n2 > 0. 

The requirement on the boundaries is that they must guarantee dynamical localisation 
(see section 121). This happens if the wave function on the boundaries is "so" small that 
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the kicking potential cannot spread a "substantial" part of the wave function out of 
the boundaries. The compatibility of the values of the parameters involved - t, k, ui, n2 
- is checked using a consequence of the conditions that grant a fractional dimension 
of the graph of the survival probability (see section |2I). This consequence is that the 
square of the wave function decreases with time keeping its shape constant, in the limited 
momentum lattice representing the open system. 

Let us recall that the typical shape of a one-dimensional localised wave function 
is exponential, extending in a region intermediate between the support of the initial 
state in momentum space and the absorbing boundary. In a linear-logarithmic plot the 
wave function is (apart from erratic fluctuations around its mean decrease) a line in this 
intermediate region; constancy of the shape means constancy of the steepness of the 
line. This criterion, which is in fact a localisation criterion, was used as a prerequisite 
for all our simulations. If the boundaries are too far away from the initial state, the 
decay is extremely slow (a consequence of strong dynamical localisation). To avoid long 
waiting times (which are hard to reach experimentally), asymmetric boundaries have 
been used, with 1 ^ |ni| <^ |n2|, and a statistical initial ensemble of orbits at t = with 
p = and randomly distributed phases 6, i.e., ipin^t = 0) = 5ri,o- The wave function 
in momentum space ip{n;t) then evolves to a shape which is asymmetric with respect 
to n = 0. On the side where the cutoff is closer to the origin, the wave function does 
not decrease exponentially, and in a linear-logarithmic plot the momentum distribution 
shows a broad and smooth maximum, while at = a sharp peak would be present 
if we choose 1 ^ |ni| ~ \n2\. Although the exponential decrease on the side of the 
larger cutoff n2 is influenced by the opening at ui, the shape indeed remains constant 
for a sufficiently large number of kicks in a range [fi, Ti2[, where the precise value of 
n ^ 50 . . . 100 depends on the choice of Ui. 

4. Numerical Results for fixed quasimomentum 

The central result of this paper is the computation and fractal analysis of the survival 
probability Psurv(T; /3, /c, ni, 712) as a function of r, while the other parameters are 
fixed for each curve. Our fractal analysis comprehends the computation of (a) the box- 
counting dimension [71 (b) a variational algorithm dimension |^, together with 
the calculation of the (c) correlations, and (d) variances of the graph Psurv(T). Several 
curves are computed with different choices of parameters. Our results are essentially 
independent of quasimomentum j3 and the applied boundaries ni and n2, whose choice 
is guided by the considerations stated in section 13.31 

Numerical algorithms, of course, do not distinguish the origin of the irregular profile 
of a fractal graph. To make sure that the observed fractality is actually produced by 
quantum effects, we verified that the increase of k in the range 2 ... 6 (for r > 1.4) 
is accompanied by a monotonic increase of the fractal dimension. This is a signature 
of fractality owing to dynamical localisation of a weakly open quantum system. As 
k reaches a certain saturation value /csat ~ 4.5 (where kr > 5 exceeds the global 
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Figure 1. The survival probability as a function of r for fc = 5, /3 = 0, rii = —1, 
n2 — 200 and different kick numbers t. The magnification in the lower panel shows 
that, as t increases, self-affinc fluctuations occur on finer and finer scales in t. 

chaos border [201 HI] and quantum chaos is fully developed) we verified that the fractal 
dimension ceases its substantial growth observed in the range k = 2 . . . 4.5. 

At fixed kick number t, the survival probability is in principle a smooth function 
of r on a sufficiently small scale St. After a grid in r is chosen, fractality is expected to 
increase as t increases due to the appearance of fiuctuations on finer and finer scales. A 
finer grid requires a longer time to yield a "fractal graph" down to finer scales, because 
it takes longer for the fiuctuations to appear on a scale smaller than the grid resolution. 
This scenario, where fractality is generated by "dynamical intrusion" , is exemplified in 
figure H where the survival probability in the localised regime is shown after various 
interaction times. The calculation of the fractal dimension as a function of time shows 
a monotonic increase from unity up to a value between 1.6 and 1.7. 

We computed the survival probability Psurv(T) for various interaction times of up 
to 10^ kicks. The latter value is much larger than the kick numbers of the order 100 
typically realised in state-of-the-art experiments jHl Oil- Nevertheless, the monotonic 
behaviour in time can itself be used as an important signature of fractality. In figure |21 
(a-c) the profile of Psnwi'T'i t) is shown along with a small, yet representative part of three 
successive magnifications over two orders of magnitude in the kicking period r. The real 
parts of the quasienergy spectrum are presented in figure |21 (d-f) in the same ranges of 
r. The visible avoided crossings are a consequence of quantum chaotic dynamics and 
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Figure 2. (a,b,c) show the survival probabihty of figure Rafter t = fO"' kicks at 
different magnifications. For the same parameters, (d,e,f) show the real parts of the 
quasienergies as a function of r (obtained as the eigenphases of the evolution operator 
1^. which was represented in the basis of momentum states as a finite matrix in 
the range n &]ni,n2[ and then diagonahsed). We see that the fluctuations on finer 
and finer scales are accompanied by ubiquitous avoided-crossings in the eigenvalue 
spectrum (note that for better visibility in (d-f) only a small part of the full spectral 
range [— tt, tt] is shown). 



their ubiquitous presence on different scales in r naturally compares to the self-affine 
ffuctuations of the survival probabihty. This comparison highlights the fact that the 
observed fractality is indeed a consequence of quantum chaos. 

The box-counting plot in figure El (a) shows the number of adjacent squares N{6) 
of width 6 along the r axis necessary to box all points of the curve from figure El 
(a). The scaling law N{6) ~ 5"^' thus determines the fractal dimension Df. The 
variational method (b) is a substantial refinement of the box-counting which typically 
gives more reliable results • It involves the division of the full analysed r interval in 
R subintervals, and the total variation of the curve on groups of 21 adjacent subintervals 
is computed. The average of these quantities is called Vr(/) and the value of R which 
gives the best scaling of the form Vr(Z) ~ is used. In addition to the direct 
fractal analysis of Psurv(T), we computed the autocorrelations and the variances of the 
fluctuating graphs. The correlations C(Ar) = (Psurv(T) ■ -Psurv(T + At))t are shown 
in figure El (c), the variances V{At) = (|Psurv(T + Ar) — Psurv(T)P)r in figure El (d). 
Recalling the power-law scaling of p(r) (see section EI), we can check the following set 
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Figure 3. Fractal analysis of the survival probability from figure El (a) using the 
following methods: (a) box-counting, (b) variational method, (c) correlations, (d) 
variances. The exponents of the fits (solid lines) are I?f = 1.6 (a) and 1.7 (b), Ocorr = 0.8 
(c) and Ovar = 0.8 (d). 



of relations: 

p(r)~r-" (8) 

and 

C(Ar) - C(0) ~ At", V{At) ~ At" with A = 2 - a / 2 , 

in the presence of the numerically confirmed identity between the temporal decay 
exponent of Psurv(^) oc t~" and the exponent of the correlations [T21I22]- These relations 
can be used as alternative and independent routes to the determination of the fractal 
dimension D{. This follows from the fractional Brownian motion nature of Psurvl^), 
which itself originates form the spectral properties of the opened 5-kicked rotor 7j, and 
which determines the At —>■ properties of statistical quantities such as correlations 
and variances [32] ■ 

Table 1 reports the fractal dimensions which were obtained by the above four 
methods. -Df be and D{^^ are the box-counting and the variational dimension, 
respectively, while the exponents of the correlations and variances are denoted Ocorr 
and Ovar- The table highlights the features already mentioned, i.e., the increase of D{ 
for increasing kicking strength k and its basic independence of quasimomentum and the 
chosen cutoffs. The fractal dimension saturates for k > /csat ~ 4.5. We verified this 
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Table 1. Fractal analysis of the survival probabilities after 10'* kicks. I?f,bc states to 
box counting dimension, while Df ^ is obtained via the variational method. Ocorr and 
Cvar are the exponents of the fits to the correlations and variances, respectively. The 
estimated uncertainty derived from our fits is ±0.1 for the fractal dimensions as well 
as the exponents. 

saturation with a series of simulations conducted for 13 values of G [2,6] (not all 
shown in table 1). 

The obtained four independent methods of our fractal analysis (summarised in table 
1) give fairly consistent results with each other, with an estimated precision of ±0.1. 
A systematical underestimation by box-counting method is observed, but also expected 
p?T] when applying it to curves with D ^ 1.5. 

For our choice of the grid in r [see Eq. ©] we noticed by inspecting the correlations 
and variances that, for /c ~ 5, not all the fluctuations of the true curve are resolved by 
our grid. This yielded systematically smaller and meaningless values for a, a problem 
which does not affect the box-counting and variational method that do not depend so 
critically upon the values of neighbouring points of the analysed graph. Augmenting the 
resolution of our grid in r on a test interval [tq, ro + 10^ 6t] (c.f. Eq. ((Tj) for the definition 
of To and St) we nevertheless were able to estimate the exponents of the correlations 
and the variances for k > 5 and n2 > 250 shown in table 1. 

As a final test of our hypothesis that no trace of quantum resonances can be 
observed for the chosen interval in r and our maximal interaction time of 10^ kicks, 
we analysed the survival probability for k = 5 for two different quasimomenta /5 ~ 1/3 
and (3 = 0.378942469767714 (stated as 0.33 and 0.38, respectively, in table 1). The latter 
value was chosen as a fraction of the golden mean to avoid any resonance condition in (5 
[see Eq. Q]. As can be seen from table 1, no dependence on quasimomentum is found 
for the dynamically localised regime {k = 5). 

In this section we presented a full-featured analysis of the fractal dimension of the 
survival probability Psurv(T), studied the dependence on the parameters t and k and 
observed how these dependences provide systematical signatures of fractality caused by 




Figure 4. (a) shows two survival probabilities for fixed f3 « 0.006554 and /3 « 
0.002009 together with an average of 10'^ (3 values equally distributed in [0, 0.01] at 
t — 500, k = 4.5, ni = —1, n2 — 200. The average curve (thick) is smoother but its 
fractional dimension is nevertheless greater than unity, (b) shows the fractal analysis 
by the variational method for f3 « 0.006554 (circles) that yields D[ 1.6, for the 
average of 10"^ /? G [0,0.01] (diamonds) with Df « 1.2, and for the average of only 10 
values of/? in the same interval (squares) that gives the intermediate value of Df « 1.4. 



quantum effects. We found tliat Psurv(T) is indeed fractal over a substantial range of 
scales, and its dimension can be estimated between 1.6 and 1.7. These numbers are 
stable when varying the initial quasimomentum (which is a constant of the motion) 
and the selected locations of the cutoffs rii and n2. Having in mind that the numerical 
determination of the fractal dimension of a graph bears some finite error, our results 
are consistent with the fractal dimension 1.5 found for fixed r = 1.4 in the scan of 
quasimomentum ^2]- The tendency towards a slightly larger fractal dimension in our 
data could be related to the distribution of decay widths, whose precise form is sensitive 
to the chosen values of rii and n2 (see Ref. |21j and discussion in section |21). 

5. Signatures of fractality for realistic experimental conditions 

5.1. Experimental control of parameters 

To realise an experiment where the fractal dimension of the survival probability, as 
studied in the preceding section, can be measured, it is necessary to address some 
principal problems of atom-optics kicked rotor experiments. 
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Figure 5. (a) average of the survival probability (for k = 4.5 and after 500 kicks) 
over 10'^ values of /? uniformly distributed in [0,0.01] (solid) and [0,0.025] (dotted) 
together with a survival probability for a fixed value /3 ~ 0.65. (b) Same as (a) for 10'^ 
values of /3 uniformly distributed in the shown intervals, (c) Fractal analysis by the 
variational method for the survival probabilities shown in (a,b) . The fractal dimensions 
are obtained by linear fits (shown only for A/3 — 1) through the symbols Df « 1.6 
(inverse pyramids, f3 « 0.65), 1.2 (circles, P & [0,0.01]), 1.2 (squares, /3 G [0,0.025]), 1.3 
(diamonds, (3 e [0, 0.05]), 1.4 (pyramids, (3 G [0, 0.1]), and 1.5 (left triangles, /3 G [0, 1]). 



Control over the kicking strength k is granted with a precision of a few percent 
im ITij . Anyway, table 1 tells us that a variation of k of the order up to 25% is 
not crucial. Time is one of the best controlled experimental parameters, and this 
feature makes it an ideal candidate for implementing an experiment to search for 
fractal fluctuations. Kicking periods between about hundred nanoseconds and a few 
hundred microseconds are available, with a maximal precision of a few nanoseconds 

El ^1 For caesium atoms, this range corresponds to dimensionless kicking 
periods (see section ^ r 10~^ . . . 18, and a maximal precision of 6t > 10"''. This 
precision implies that about 100 points could be scanned in our analysed interval in r, 
which would be sufficient for a rough, quahtative verification of our predictions. 

Any experiment will have a finite detection window of observable momentum 
classes. The actual width of this window is typically determined by the imaging 
resolution and by the minimal signal-to- noise ratio of the measurement device |m il4[ l7rj. 
The detection window also determines a maximum interaction time after which the 
detection of a constantly decreasing atomic ensemble (due to the open boundary 
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Figure 6. (a,c) average of the survival probability over 10'^ values of (3 uniformly 
distributed in [0,0.1] and [0,1], after (a) 100 and (c) 200 kicks and for k = 4.5. The 
enhancement of self-affinc fluctuations with time is clearly visible, (b) and (d) show the 
fractal analysis by the variational method for (a) and (b), respectively, corresponding to 
/3 ~ 0.65 (inverse pyramids), or 10^ values of f3 uniformly distributed in [0.01] (circles), 
[0.025] (squares), [0,0.05] (diamonds), [0,0.1] (pyramids), [0,1] (left triangles). The 
fractal dimensions are Df « 1.5, 1.2, 1.1, 1.2, 1.2, 1.4 in (b) and 1.6, 1.2, 1.1, 1.2, 1.3, 
1.5 in (d) for increasing width of the /3 distribution. 

conditions) becomes meaningless. In other words, the maximum number of kicks 
is limited by the precision disposable in the determination of the final momentum 
distribution. Correspondingly, in our results reported below we choose the minimal 
kicking strength k = 4.5 where the fractal dimension starts to saturate (see section 
and a maximum interaction time of 500 kicks. The latter implies that we can choose a 
wider grid in r because very fine structures do not develop for interaction times t < 500. 
We used St' = 10 St, Tq = tq and m' = m / 10 [c.f. Eq. A problem will certainly 
be the realisation of our idealised absorbing boundary conditions at specific momentum 
classes of the atoms. Here methods using, for instance, external cutting potentials - 
such as so called radio-frequency knives [23] or equally operating additional lasers - 
could be thought of. 

5.2. The experimental initial ensemble 

To approach real experimental scenarios, we shall analyse the survival probability for a 
smaller number of kicks of order 100 [HI El] and take into account an initial spread of 
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quasimomentum among the ensemble of cold atoms fUUHZj. 

For a typical ensemble of cold atoms, the momentum distribution is Gaussian-like, 
with a width exceeding that of the Brillouin zone 2 hkL, equal to 1 in our dimensionless 
units Cni UH UH I2H1 ■ Folding produces approximately a uniform distribution 
in the entire Brillouin zone, i.e., a uniform distribution of quasimomenta with a width 
of A/9 = 1 pni- Using atoms in the Bose-Einstein condensate phase as initial ensemble 
allows the experimentalist a much better control over the width of the quasimomentum 
distribution jHSI- Values of A/? ^ 0.05 have been realised in this context [IHl EHl 1^ . 
Letting the condensate expand a little before the actual kicking evolution, allows one 
to reduce the atom-atom interactions to negligible values, with only slight changes 
in A/9 As a consequence, the survival probability, experimentally measured by 

counting the number of atoms contained within the finite detection window, would be 
the result of an average of many independent survival probabilities with different values 
of quasimomentum. The independence of probabilities follows from the independent 
dynamics of the atoms [HI CHI HH CH EZl IHE] , while the coherent evolution of a single 
atoms is still essential for the observed behaviour. 

We computed Psurv(T; P) for different ranges of /3, and then averaged the resulting 
curves to arrive at (-Psurv(T; /?))/3- Figure El investigates the effect of averaging over /3 
on the fractal dimension. The survival probabilities for two fixed (3 are shown, together 
with the average for a uniform distribution of 10^ values of /5 G [0,0.01]. Figure ID 
(a) shows that the average curve is quite smooth on large scales, but nevertheless 
presents fluctuations on finer resolutions, with a fractal dimension substantially larger 
than unity. We verified that, by decreasing the number of atoms in the ensemble, the 
dimension steadily increases. We encounter a signature of fractality, which experiments 
could detect even far from the idealised limit of the one-atom dynamics. That is, the 
fluctuating behaviours of the averaged curves is a direct consequence of fluctuations of 
single P curves. 

Figure El (a) and (b) show the average survival probability for ensembles with 
the same number of (3 values but with different, larger widths A/3 of the initial 
quasimomentum distribution. Wider distributions are smoother on large scales and are 
not drawn in figureEl(a) because they could not be appreciated by eye when compared to 
curves for A(3 = 0.010 and 0.025. The magnification in (b) shows that the fluctuations 
exhibit smaller excursions. A fractal analysis [see figureEKc)] by the variational method 
shows that the dimension Df remains in all cases larger than unity and, moreover, does 
not vary monotonically as AjS is increased. 

We interprete our results for finite A/3 in the following way: while a small range 
Af3 tends to wash out the fractal behaviour of the curves with one fixed /?, an average 
over larger ranges Af3 tends to lift the fractal dimension again. This line follows nicely 
from the prediction of [Z] where it is argued that fractality can arise from superimposing 
non- fractal patterns on appropriate scales of the scanned variable [221123 • 

Figure El repeats the analysis of figure El for t = 100 (a,b) and t = 200 (c,d). The 
fractal dimension of each average survival probability for a definite value of A/3 is seen to 
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Figure 7. (color online). On the left panel the survival probability as a function of 
/? and T is shown after 500 kicks, for k = 4.5, n\ ~ ~\ and 712 = 200. Slices of the 
graph Psurv(/3,T) are shown on the right panel. The thick lines represent the survival 
probability as a function of r that is analysed in section 0] (we have used 10'^ similar 
curves to compute the incoherent average which is the experimental observable, as 
explained in section [5.2(1 . The thin curve lying in the plane orthogonal to the t axis is 
the survival probability Psm v (/3) as a function of the quasimomcntum /3 such as studied 

in [ig. 

be a monotonic function of time, what points out once more the dynamical origin of the 
analysed fluctuations. This contrasts the dependence of the fractal dimension on A/? at 
fixed time, which is non monotonic, because averaging both washes out the fluctuations 
of the single curves for small A/5 ~ 0.1, while it creates new ones by superimposition 
for A/? > 0.1. 

From a more general perspective, the survival probability Psurv(/5, ''") can be seen as 
a surface lying over the plane spanned by the two variables r and /3. In ^21 -fsurv(/?, t) 
was analysed at fixed r using (5 as scanning parameter, i.e., a slice of the surface parallel 
to the (5 axis was analysed. In section 0] we studied the "orthogonal" problem of fixed 
/3, using r as a scanning variable. Averaging over /? can be interpreted as a "column 
density", i.e., Psurv(/?, t) integrated over the /? degree of freedom. The resulting averaged 
curve is the observable experimentally accessible as discussed in this section. Looking 
at figure our results can thus be interpreted geometrically: the fractal behaviour of 
the slices Psurv(''"; /5 = const) is to a large extent preserved by the average over a typical 
experimental spread in /3. 

6. Conclusions 

We considered the quantum kicked rotor, a paradigmatic model of quantum chaos, 
which describes the time evolution of noninteracting cold atoms in periodically flashed 
optical lattices. Imposing absorbing boundary conditions allows one to probe the 
transport properties of the system, and in particular to define the survival probability 
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of atoms on a finite region in momentum space. For fixed kick numbers, tlie quantum 
survival probability depends sensitively on the parameters of the system, and a self- 
afiine structure of the survival probability Psurv is predicted, as either the kicking period 
or quasimomentum is scanned. 

Instead of using the initial quasimomentum (3 as control parameter, as done in the 
numerical simulations of ^21; we used the kicking period r as scan parameter, which 
is much better controllable experimentally. We verified the fractal nature of the graph 
of the survival probability Psurv (t) in the dynamically localised regime, and obtained a 
fractal dimension D{ ^ 1.6±0.1 for large but finite interaction times, for which quantum 
resonances do not manifest. 

Any experimental setup prepares cold atoms with a finite spread in quasimomen- 
tum. The experimental observable is then the average of the survival probabilities over 
the quasimomentum distribution. We reproduced this observable by computing the 
incoherent average (Psurv (t))/?, and found that the fractal dimension of the average re- 
mains substantially larger than unity, even for shorter interaction times of a few hundred 
kicks. 

We conclude that the fractality in the survival probability induced by quantum 
chaos is an unexpectedly robust feature and in spite of many challenging aspects (see 
section could be observed in a future atom-optics experiment. Apart from the 
experimental verification of fractal fluctuations of purely quantum origin, a remaining 
open problem is whether a universal scaling law for the fractal dimension could be found 
as a function of both parameters r and f3, including quantitative predictions for the here 
computed averages (Psurv (t, /?))/? over a finite range of j3. 
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